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It is now widely accepted that most of mass— energy in the universe is unobserved except by its gravitational effects. 
Baryons make only about 4% of the total, with "dark matter" making up about 23% and the "dark energy" responsible 
for the accelerated expansion of the universe making up the remainder. We focus on the dark matter, which is the 
primary constituent of galaxies. We outline the observed properties of this material, enumerating some candidates 
covering 90 orders of magnitude in mass. Finally, we argue that the weak scale (100 GeV) is relevant to new physics, 
including the dark matter problem. 



1. THE DARK MATTER PROBLEM 



O 
O 

<N 

Cosmological parameter estimation, based on cosmic microwave background anisotropies measured by the WMAP 
' satellite 0] and based on the power spectrum of galaxy density fluctuations measured by the SDSS collaboration |3] 
indicates that most of matter in the universe is unobserved save by gravity. Approximately 4% of the cosmological 
energy density is accounted for by baryons, 23% by the dark matter discussed in this lecture, with the remainder 
being the "dark energy" responsible for the observed accelerated expansion of the universe. The dark matter is 
the primary component of cosmological structures at the dwarf galaxy scale and larger. The dark matter and dark 
energy are of unknown composition. In this lecture we discuss the physics of dark matter, and enumerate a list of 
, candidates, not meant to be exhaustive. 

o 



2. DARK MATTER PROPERTIES 

In this section we enumerate five basic properties that dark matter is observed to have. The first three (that it must 
be non-radiating, approximately collisionless and non-relativistic) do not otherwise place any positive constraints on 
the space of possibilities, while the last two (that it must be fluid and classical) place upper and lower bounds 
C$ ■ respectively on the mass of the particles. 

> : 

X ' 2.1. Optically Dark (Dissipationless) 



Dark matter is not observed to shine, thus the dark matter particles must have very weak electromagnetic in- 
teractions. Either the electric charge and electric and magnetic dipole moments must vanish (or be very small), 
or the particles must be very heavy. For a review of constraints on dark matter from background light across the 
electromagnetic spectrum, see Ref. An important consequence of this is that the dark matter can not cool by 
radiating photons, and thus will not collapse to the center of galaxies as the baryons do, by radiating their energy 
away electromagnetically. In other words, the dark matter is very nearly dissipationless. 



2.2. Collisionless 

In addition to not interacting with "light" matter (to ensure that energy is not dissipated), the dark matter must 
be nearly collisionless as well. Even if the dark matter can not radiate energy, collisions will serve to make halos 
round, in contrast to data that overwhelmingly indicates triaxiality, e.g. in clusters The limit of this constraint is 
that there can be dark matter self interactions that are important at high densities and short distance scales, which 
may serve to erase small scale structure in galaxy halos, and erase the density cusps expected at the cores of galaxies. 
These effects appear for interaction cross sections of a ~ (m/GeV)(A/Mpc)~ 1 barns U where A is the mean free 
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path of the particles in typical halos. Essentially, this order of magnitude is indicative of hadrons, and there are a 
number of possibilities. 



2.3. Cold 

In the early 1980's it was noted that "cold" dark matter explained the observed properties of galaxies quite well 
0. The measured two point correlation function of galaxies indicates that there is a large amount of power on 
small scales. If the dark matter particles have significant velocities, then the small scale power is erased. For this 
not to occur, the particles must be sufficiently non-relativistic at the epoch of matter-radiation equality, when the 
temperature of the universe was roughly 1 eV Q- In practice, the constraint is that a particle species in thermal 
equilibrium must have a mass larger than about 1 keV [a|. Of course, a non-thermal species can have a smaller mass. 
Dark matter that is not utterly cold may in fact ameliorate some of the difficulties of the cold dark matter model 
:9j; some suppression of small scale density fluctuations may be favored. 



2.4. Fluid 

The dark matter must be sufficiently fine on galaxy scales such that the discreteness is not detectable as yet. 
There are two basic effects to be concerned with. The first is that the granularity of dark matter provides a time 
dependent gravitational potential, which may disrupt bound systems. For particles of 10 6 Mq, this effect would heat 
the galactic disk to observable levels, while globular clusters would be disrupted for smaller masses, around 10 3 Mq 
Second, the discreteness of dark matter introduces Poisson noise in the power spectrum of density fluctuations, 
which conflicts with observations of the Lya forest if the particles are more massive than roughly m > 10 4 Mq [llj. 
We thus find that an upper bound on the mass of particles of roughly m < 10 3 ~ 4 Mq ~ 10 70 - 71 eV. 

Gravitational microlensing places further constraints on dark matter in the solar mass range. Searches for lenses 
by the MACHO and EROS collaborations have failed to find enough microlensing events to explain the full dark 
halo of the Milky Way [12J, though MACHO reports that 20% of the Milky Way halo consists of objects of roughly 
0.3 Mq. The constraints are that objects between 1O _7 M0 and 10 Mq do not make up more than about 20% of the 
dark matter in the Milky Way halo. Microlensing surveys are not able to place constraints much outside this mass 
range due to the limitations of microlensing event timescales. 



2.5. Classical 

Dark matter must behave sufficiently classically to be confined on galaxy scales. Bounds can be placed on the 
masses of both bosons and fermions, based only on the observed properties of galaxies: namely galaxy densities must 
reach of order GeV cm~ 3 , their velocity dispersions are of order 100 km s -1 , and their sizes are of order kpc. 

2.5.1 . Bosons 

If dark matter consists of bosons, their quantum nature is manifest only if their mass is exceedingly small. In 
order to form galaxies, dark matter particles must be confined on kpc scales. With typical galactic velocities, the de 
Broglie wavelength of the particles is A ~ (eV/ m) mm. Setting this wavelength to 1 kpc, we find that particles with 
mass m ~ 10~ 22 eV are barely confined. This mass denotes the lower end of the possible range for dark matter. In 
fact it was proposed as a mechanism to erase the small scale power in the galaxy power spectrum , called "fuzzy" 
cold dark matter. This also disallows the central density cusp found in numerical simulations of galactic dark halos. 

2.5.2. Fermions 

The lower bound on the mass of a fermionic dark matter particle is much more stringent than that for bosons. 
As first pointed out in Ref. [14|. the phase space density of fermions has a maximum value / = gh~ 3 , where g is 
the number of internal degrees of freedom, and in fact for a relativistic gas of fermions in thermal equilibrium it is 
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half this value. In a galactic halo, the maximum phase space density assuming a Maxwellian velocity distribution is 
simply / = p/[m 4 (27T(T 2 ) 3//2 ]. Pauli blocking then enforces m 4 > ph 3 /[g(2na 2 ) 3 ^ 2 ]. Taking a Milky Way type galaxy, 
where p > 1 GeV cm~ 3 is required in the center and a = 150 km s _1 , and assuming g — 2, we find m > 25 eV for 
fermions. 



The dark matter problem is a significant challenge for the particle physics community in that any candidate must 
satisfy stringent constraints. From a particle physics point of view, the properties discussed in section|21are not hard 
to arrange, but two implicit properties are: namely that the candidate have a lifetime much longer than the Hubble 
time ~ 10 Gyr, and also that its cosmological density is compatible with the observed dark matter density. The 
current measured value of the dark matter density is JIcdm^ 2 = 0.1351q'qq9 fljj], where ft is the density relative to 
the critical density p c = 3H 2 /(8itG), and h (= OJlj^o:?) * s tne sca l e d Hubble constant H = 100 h km s _1 Mpc -1 . 
In particle physics units, p c 10.5 h 2 keV cm -3 . Accounting for the possibility of a significant neutrino density (see 
section FO)) . we take a 95% confidence region 0.087 < ^cdm/i 2 < 0.129. 

The simplest possibility would be for a candidate to have the correct density from thermal processes alone. Com- 
puting the relic density proceeds as follows ^| . The species is assumed to be in thermal equilibrium at early enough 
times, i.e. when the universe was essentially a radiation dominated plasma. The production and annihilation rates 
are tracked as the universe expands, and when these rates become longer than the (current) Hubble time, the species 
is said to be "frozen out" with a comoving density that no longer changes. This statement is expressed in the 
Boltzmann equation, 



where n and n cq are the particle density and its equilibrium value, H is the Hubble parameter and (av) is annihilation 
cross section, appropriately thermally averaged [16j. We immediately notice that in the absence of interactions, the 
number density evolves exactly as expected, n ex a~ 3 , where a is the cosmic scale factor and H = a/ a. The relativistic 
gas has a temperature evolving as T oc a -1 , thus n cx T 3 . The details of freeze out now depend on whether the relic 
is relativistic or not. 

3.1. Relativistic Thermal Freeze-out 

For relativistic particles, n cq oc T 3 , thus a species that has frozen out will still track the equilibrium density. The 
relic density of such particles is thus insensitive to the details of freeze-out. The only relevant parameter is the 
effective number of degrees of freedom at freeze-out, g+. This can change as a function of temperature, thus the 
density of a frozen out species need not track its equilibrium density, but the scaling with temperature remains as 
long as the particles are relativistic. The well known example of neutrinos freeze out with g+ ss 10.75, prior to the 
freeze-out of electrons and positrons (which lowers <?*). Thus the (massless) neutrino temperature is lower than the 
photon temperature in the present universe. The relic density of any relativistic species is given by 



3.2. Non-Relativistic Thermal Freeze-out 

Non relativistic particles have an exponentially suppressed equilibrium density, n cq oc (mT) 3 / 2 exp(— m/T). A 
frozen out species would thus have a density much larger than its equilibrium density. We can simply derive the scaling 
of relic density with annihilation cross section as follows. In a radiation dominated universe, H oc T 2 (since H 2 oc p 
according to the Friedman equation). The freeze-out condition can then be expressed simply as H(Tf r ) ~ ero oq (Tf r ). 



3. COSMOLOGICAL RELIC DENSITY OF PARTICLES 




(i) 




(2) 
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We find that Tf r oc m, with only logarithmic corrections (this is due to the exponential Boltzmann factor). We can 
then derive the current mass density as p ~ mn cq (Tf T )(T/T{ v ) 3 ~ T 3 (m/Tf r )/cr. Thus we find the crucial result 
that p oc 1/cx, with only logarithmic dependence on mass since m/Tf r is roughly constant. This means that smaller 
annihilation cross sections yield larger relic densities, which makes sense, as less efficient annihilations should allow 
more particles to survive. Taking the simple scaling that a oc m~ 2 , we find the approximate relic density 



nh ~ liwJ ■ ( 3 ) 

This is a very interesting result in that a stable particle at the weak interaction scale of several hundred GeV would 
give the proper relic density to be dark matter. Any Weakly Interacting Massive Particle (WIMP) might thus be a 
compelling dark matter candidate. 



3.3. Pedagogical Example 

We can illustrate both the relativistic and non-relativistic freeze-out regimes with a single toy model as follows. 
Taking a Dirac fermion annihilating through a wide Z' boson, the annihilation cross section is 

9 4 



aV \Mvcos A O w ) (s-m%,) 2 + m%,' (4) 

In figure ^ we plot the relic density of this particle for a wide range of masses, as calculated by the DarkSUSY 
code |17|. Three regimes are evident. For masses below 1 MeV, the freeze out is relativistic, and the relic density 
is proportional to mass. Above 1 MeV, but below mz 1 the interaction is like low energy weak interactions, with 
a oc E 2 = m 2 , thus the relic density is proportional to the inverse square of the mass. For masses larger than the Z', 
the usual behavior of a oc m~ 2 is recovered, with relic density proportional to the square of the mass. In a parallel set 
of three regimes, the boundaries between hot, warm, and cold dark matter are approximately illustrated. Comparing 
with the known value of the relic density, the hot and cold possibilities appear at 10 eV and 1 TeV respectively. 




og (m x /GeV) 



Figure 1: Relic density of a species freezing out in thermal 
equilibrium. The interaction cross section is of weak interac- 
tion strength. For masses below ~ MeV, the particles are rela- 
tivistic at freeze-out, above this mass they are non-relativistic. 
Above a mass of ~ 100 GeV, the interactions look "electro- 
magnetic" as the gauge boson becomes light relative to the 
particle mass. The three regimes of cold, warm, and hot dark 
matter are illustrated, along with the WMAP 2er constraint 
on dark matter density. 
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4. DARK MATTER CANDIDATES 



In this section we outline some of the proposals for dark matter. While this list is by no means exhaustive, we 
will attempt to cover the range of possibilities that have been considered at least qualitatively. We will proceed 
very roughly from the smallest mass candidates at 10 -22 eV (the fuzzy cold dark matter discussed previously) to the 
largest at 1O 6 M « 10 72 eV. This range covers more than 90 orders of magnitude in mass; unfortunately observations 
do not place stronger bounds on possible candidates. However, the weak scale (100 GeV) is quite interesting for the 
dark matter problem, as we will describe. 

4.1. "Fuzzy" Cold Dark Matter 

This model was discussed in section I2. 5. II 

4.2. Chaplygin Gas 

The Chaplygin gas is an attempt to unify dark matter and dark energy in a single fluid with an equation of state 
p = — A/ p. Generalizations of this have also been considered This type of fluid can arise in certain string- inspired 
models involving d— branes. The transition between the dust dike universe at moderate redshift to the accelerating 
universe today occurs naturally in these models. 

4.3. Axions 

Strong interactions naively violate CP with order unity strength (the coefficient 9 of the FF term in the QCD 
Lagrangian is not obviously small), while experimentally the violation must be smaller than 8 ~ 10~ 9 . Promotin g th e 
parameter to a dynamical field solves this problem, at the expense of a new light pseudoscalar boson, the axion [19J. 
The cosmology and astrophysics of these particles has been studied at length, beginning with Ref. |20| and reviewed 
in Refs. |2lj . To summarize, there is an allowed mass window of interest to dark matter, with fieV < m < meV. The 
lower bound comes the fact that the relic density of axions Q a h 2 ~ (4/ieV/m a ) 7 / 6 (200MeV/A QCD ) 3/4 , while the 
upper bound comes from the fact that the axion interactions must be weak enough that stars cool via photons and 
not axions. Axion production is non-thermal, thus they are cold, even though their masses are far below the ~ keV 
limit for warm dark matter. Early on it was pointed out that axions could be detected by resonant axion-photon 
conversion in a magnetic field [22j | , an experimental program that has been undertaken with some success |'2^ . Axion 
dark matter is particularly interesting from an experimental point of view in that the velocity distribution of the local 
dark can be measured with high accuracy. A record of the process of virialization would be evident, of considerable 
interest to studies of galaxy formation. 

4.4. Neutrinos 

We now know that neutrinos are one component of the cosmological dark matter. Neutrino oscillation experiments, 
both for atmospheric neutrinos ^24] and for solar neutrinos [25f . have now conclusively shown that at least two of 
the neutrinos in the Standard Model are massive. Measurements of the endpoint spectrum of tritium beta decay 
place a 95% confidence upper bound on the electron neutrino mass, m Ve < 2.5 eV |2fij| . thus bounding all three 
neutrino species as the mass splittings measured by the oscillation experiments are both much smaller than this. 
The cosmological density of light neutrinos (m v < MeV) is known to be VL y h 2 = J2 "W(94.0 eV) ,27]. The WMAP 
bound 1] on this quantity is Q„h 2 < 0.0076 — > J2 m v ~ 0.7 eV, thus the cosmological bound is more stringent 
than the laboratory one (^ m„ < 7.5 eV). The measured oscillation parameters enforce that ^Z m v ~ 50 meV, and 
thus 0.0005 < fl v h 2 < 0.0076, a known constituent of dark matter as even the bottom of this mass range would be 
non-relativistic today. However, no fermion in this mass range can be a significant constituent of galaxies due simply 
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to the Pauli exclusion principle discussed above. In addition to not being abundant enough, these Standard Model 
neutrinos are "hot" and thus unsuitable to explain the properties of galaxies. 

Very massive neutrinos have also been considered as dark matter candidates, first in Ref. 28]. We now know that 
none of the Standard Model neutrinos are suitable, but their right-handed counterparts might be acceptable. Active 
neutrinos in a 4th generation would be acceptable except for the fact that they would have been detected in sensitive 
underground experiments by now. 

4.5. MeV Dark Matter 

The INTEGRAL satellite has observed an excess of 511 keV positron annihilation emission toward the galactic 
center, which is difficult to reconcile with astrophysical sources. It has been suggested that this emission might be 
due to e + e~ pairs that are the annihilation products of dark matter particles in the MeV range j2Q|. Their mass 
must be low to exclude hadronic final states, which would overproduce photons from e.g. tt° decays, and also so that 
the leptons are not so energetic that they escape the galactic center region. However, the particle physics motivation 
for this scenario is weaker than for other candidates. 



4.6. Supersymmetry 

Our current understanding of physics includes no fundamental relationship between bosons and fermions. Super- 
symmetry was proposed as just such a relationship, first on the 2-d worldsheet for string theory |3lJ, and slightly 
later for standard 4-d spacetime Q Supersymmetric theories put bosons and fermions in common multiplets; for 
most such theories then the corresponding states must have similar masses. It is thus obvious that if supersym- 
metry is a valid symmetry, it is badly broken in the low energy world in which we live. These theories admit the 
so-called non-renormalization theorem, which remains valid if the supersymmetry is broken "softly" . Very roughly 
stated, the only infinite renormalizations are of wavefunction type, and those are only logarithmically divergent. The 
quadratically divergent corrections in the Standard Model are absent. 

The Minimal Supersymmetric Standard Model (MSSM, for early reviews see Refs. |33j) was proposed as a model 
with broken supersymmetry. Each Standard Model state has a superpartner with spin differing by h/2: matter 
fermions (quarks, leptons) have scalar partners, gauge bosons have spin 1/2 partners, and Higgs bosons (two SU{2) 
doublets are required) also have spin 1/2 partners. The soft supersymmetry breaking in this model allows for 
the masses of superpartners to differ from their Standard Model counterparts. The most general gauge invariant 
Lagrangian involving these fields includes terms that violate either baryon or lepton number, e.g. a term allowing 
scalar up-type squarks to decay directly to a down quark and charged lepton. If both baryon and lepton number 
violating interactions are included, proton decay occurs at the weak scale. The usual solution is to impose R-parity 
conservation, where superpartners have negative R-parity and thus must interact in pairs. This means that the 
lightest superpartner (LSP) is absolutely stable. Early on, it was realized that the LSP (first considered to be the 
gravitino) was thus a dark matter candidate [34j . 

4.6.1 . Gravitinos 

The superpartner of of the graviton, the spin 3/2 gravitino, was the first SUSY particle considered for the dark 
matter problem [34j |. In models where the gravitino is the LSP, it is often quite light (keV), and would thus be warm 
dark matter. In cosmology, the overproduction of gravitinos is somewhat problematic, though not insurmountably 
so Is^. Gravitinos at the weak scale, whose relic density would be obtained through the decays of the next lightest 
superpartner, are also an interesting possibility puj . 

4.6.2. Neutralinos 

The favored supersymmetric dark matter candidate is the lightest neutralino; these are the four spin 1/2 Majorana 
fermion superpartners of the neutral gauge and Higgs bosons (usually denoted X1-4) |3?j | . Similarly, there are two 
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charged Dirac fermion superpartners of charged gauge and Higgs bosons, the charginos xt-2- The successes of 
supersymmetric models depend on the fact that the gauginos have weak-scale masses, thus the proper relic density 
for any stable states comes essentially for free. We focus on neutralinos as dark matter in section 

4.6.3. Sneutrinos 

The scalar partners of the neutrinos are possible dark matter candidates, but they are disfavored for two reasons. 
First, they annihilate quite efficiently, requiring masses above 500 GeV to provide a relic density consistent with dark 
matter. Second, their elastic scattering cross sections on nuclei are quite large (on the order of femtobarns), to the 
point that they would be easily detectable in current experiments |38( . 

4.6.4. Axinos 

If the axions proposed to solve the strong CP problem exist (see section l4~3f and in addition supersymmetry is 
valid, the axion will naturally have a spin 1/2 partner, the axino. Depending on the conditions in the early universe, 
these might be either warm or cold dark matter |39| . 

4.6.5. Q-balls 

Supersymmetric theories typically permit non-topological solitons, dubbed Q-balls. These carry baryon and/or 
lepton number, and can be absolutely stable if large enough. They are an interesting dark matter candidate, and 
may even have relatively strong self- interactions at the level discussed in section 12721 40] . 

4.6.6. Split SUSY 

Recently it has been noted that the successes of supersymmetric models hinge on the fact that the gauginos 
have masses at the weak scale, while having the scalar superpartners at the weak scale is somewhat problematic. 
If the Higgs hne-tuning problem is simply ignored (as the cosmological constant fine-tuning often is), the scalar 
superpartners can be made very heavy, keeping light gauginos and Higgsinos [4lj |. This scenario has an interesting 
phenomenology, both for collider signatures and for dark matter [42| . 



4.7. Universal Extra Dimensions 

If our four-dimensional spacetime is embedded in a higher dimensional space, excitations of Standard Model states 
along the orthogonal dimensions (called Kaluza-Klein excitations) may be viable dark matter candidates Q|. For 
example, the first excitation of the B boson (associated with the U(l)y hypercharge gauge group of the Standard 
Model) has been considered. This is similar to supersymmetry in that known particles have partners, but the partners' 
spin does not differ. This scenario has been dubbed "bosonic supersymmetry" . The stability of the lightest of the 
KK excitations can be arranged by a parity symmetry, and masses around 1 TeV provide reasonable relic densities 
from thermal freeze-out, the same as with supersymmetric models. 



4.8. Branons 

String theory naturally contains objects of many different dimensions, called branes. These would naturally have 
fluctuations characterizable as particles, the so-called branons. These fluctuations can be made into suitable cold 
dark matter candidates, both thermal and non-thermal [Sf. 
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4.9. Mirror Matter 



The (modern) concept of a mirror world is an old one, extending back to the non-conservation of parity in weak 
interactions. In a mirror scenario, the dark matter could just be ordinary matter in the mirror world, where the only 
communication is gravitational. This scenario can be constructed in a braneworld context, where our universe and 
a mirror universe are two branes in a higher dimensional space. For more complete discussions of this scenario, see 
e.g. 0- 

4.10. WIMPzillas 

At the end of inflation, gravitational interactions alone can produce copious particles. For mass scales of 10 13 GeV, 
these particles, if stable, could even account for the dark matter 0|. In addition, such particles might decay with a 
lifetime much longer than the age of the universe, providing a source of ultra high energy cosmic rays - this is the 
so-called "top-down" scenario for UHECR production. 

4.11. Primordial Black Holes 

Under the right conditions, primordial black holes can form in the early universe, e.g. see Ref. 0]. Production is 
enhanced during periods where the equation of state softens (p < p/3), such as during a first order phase transition. 
This is easy to understand: if the pressure support lessens, objects collapse more easily. The last such phase transition 
in the universe is the quark-hadron phase transition, at a temperature T ~ 100 MeV. The mass contained in the 
horizon at this epoch is very roughly a solar mass. 

5. FOCUS ON NEUTRALINOS 

The lightest neutralino is the dark matter candidate that has endured the most scrutiny. In this section we 
elaborate on the discussion of the MSSM, and discuss the prospects for discovering neutralinos astrophysically. The 
MSSM starts with the Standard Model field content, adds one additional SU(2) doublet of Higgs bosons (with the 
opposite U(l)y charge), and then adds the superpartners of each of these states. The Lagrangian includes all SUSY 
conserving terms that respect R-parity, and SUSY-breaking soft terms. More than 100 new parameters are required. 
One parameter of the Standard Model is actually removed, the quartic self-coupling of the Higgs boson. These 
couplings appear in the MSSM, but with the quartic coupling Asm replaced by terms in the gauge couplings g 2 and 
g . Clearly, a theory with 100+ free parameters is not very predictive. To make some connection to reality, a simpler 
model is required. 

5.1. Minimal Supergravity 

The best studied model of low-energy supersymmetry uses gravitational effects to break supersymmetry, called 
Minimal Supergravity (mSUGRA) 0|. The model has four parameters, plus one sign, beyond the Standard Model 
ones (and again, Asm is absent). These continuous parameters are mo, the universal scalar mass at the GUT scale, 
m.1/2, the universal GUT-scale gaugino mass, A$, the GUT-scale trilinear coupling, and tan/3, the weak-scale ratio of 
the vacuum expectation values of the two scalar neutral Higgs fields. The Higgsino mass parameter |/i| 2 is determined 
by enforced radiative electroweak symmetry breaking, but the phase of /i is undetermined: it is chosen to be ±1 to 
conserve CP. 

In studying this model, a set of parameters is evolved according to the renormalization group equations (RGEs), 
checking that the gauge couplings unify, and electroweak symmetry breaking occurs. There are a number of publicly 
available software packages that solve RGEs: we use ISAJET |49|. Further tests against accelerator data (e.g. 
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sparticle masses, BR(fo — > sj), etc.) are performed. Finally, parameter sets where the LSP is not a neutralino are 
discarded. 

We can roughly understand the weak-scale parameters as follows. As the gaugino masses are unified at the GUT 
scale, their weak-scale masses are related to the ratios of the gauge coupling constants at the weak scale: strong, 
weak, hypercharge gaugino masses have M3 : M2 ■ Mi ~ 6 : 2 : 1, namely the B-ino is the lightest gaugino. Gauge 
couplings run scalar masses higher, while Yukawa couplings run scalar masses lower: squarks are heavy, sleptons are 
light, with Ul usually heaviest and tr usually lightest. 

To connect these models to cosmology, we must calculate the relic density of neutralinos. We use the publicly 
available DarkSUSY code for this purpose [r?j |. The code actually calculates the relic density of all neutralinos, 
charginos and sfermions; if the mass splittings are small, the "coannihilation" effects can be crucial. Given the current 
experimental limits, three regions of mSUGRA parameter space stand out as giving a relic density compatible with 
cosmology: stau coannihilations, where the Xi ~ tr splitting is small, the focus point region where the Xi ~ X2 ~ X\ 
splittings are small, and the region (at large tan/3) where the annihilation through the pseudoscalar Higgs A is 
resonant, XX -> -4° -> //• 



5.2. MCMC Scanning 

The regions of mSUGRA parameter space compatible with cosmology are quite narrow, and difficult to find in 
general. Grid searches of the parameter space are ~ 1% efficient at finding models passing the 95% confidence cut 
mentioned in sectional However, we can use the relic density of models to guide the search through parameter space, 
using the so-called Markov Chain Monte Carlo (MCMC) [sof . The basic idea is to use a likelihood function C = e~ x I 2 
(here taken as the WMAP likelihood function for relic density) calculated at the current point in parameter space. 
The likelihood is calculated at a proposed point, and the "chain" advances if the likelihood is higher. With lower 
likelihood, the chain advances with probability C new / Cold- In this way chains of models can be constructed that 
efficiently sample the interesting regions. A 2 dimensional example is illustrated in figure[5| In fact we have explored 
the full 4 dimensional mSUGRA parameter space, finding efficiencies with MCMC of 20%-25% (a huge improvement). 
We scanned 2.4 million models in total, finding 500 thousand passing the WMAP cut. 
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Figure 2: MCMC scan in 2 dimensions. We have fixed Ao = 0, 
tan/3 = 30, mt = 175 GeV. Black points illustrate models 
with relic density too high, red points pass the 95% confi- 
dence WMAP cut, and blue points fall below (thus they are 
acceptable, but not as a full solution for the dark matter prob- 
lem. The stau coannihilation region clearly appears at the left 
edge, while the focus point region (chargino coannihilation) is 
clear at the right edge. At larger tan/3, the Higgs resonance 
region would appear. The stau coannihilation region ends at 
m i/2 ~ 1100 GeV, while the focus point extends at least to 
mi/2 ~ 20 TeV (10 TeV neutralinos). 
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5.3. Direct Detection 



On of the most promising experimental techniques for discovering neutralinos is "direct detection" . The idea is 
that the neutralino will scatter on a nucleus, depositing energy that can be detected in a number of ways [si|. The 
scatterings are very rare (the current limits are of the order of a ~ 1CP 42 cm 2 ), and they deposit very little energy 
(typically ~ 10 keV), so the experimental challenge is considerable. In figure |3| we plot all of the mSUGRA models 
consistent with solving the dark matter problem. The axes are neutralino mass and spin-independent neutralino- 
nucleon elastic scattering cross section. The spin independent scattering amplitudes essentially couple to the mass 
of the nucleon, and add in quadrature due to coherent effects (the 10 keV exchanged is long wavelength compared 
to a nucleus). Current experimental efforts are just starting to probe the mSUGRA region, with several orders of 
magnitude improvement being quite attainable in the near future. Spin-dependent interactions are also present, but 
they are experimentally more difficult to explore, requiring odd-A targets (though the cross sections can be quite 
large). In figure 0] we plot the spin-dependent cross sections on protons and neutrons. 
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Figure 3: Direct detection of neutralinos. The spin- 
independent elastic scattering cross section of neutralinos on 
protons is plotted against neutralino mass. Red crosses indi- 
cate the mSUGRA parameter space, while green circles indi- 
cate those models that could make a significant contribution 
to the g — 2 of the muon, which may be experimentally favored 
0. Current limits from the EDELWEISS and CDMS experi- 
ments are plotted [s^l . along with the projected sensitivities of 
several proposed experiments (j^l . The region consistent with 
the DAM A annual modulation claim is illustrated as well |5fij ; 
this region is inconsistent with the negative results of other 
experiments. 



5.4. Indirect Detection 



Here, we will briefly mention the indirect techniques for neutralino searches. The first relies on the fact that 
massive objects such as the Earth or Sun will collect neutralinos due to elastic scattering into bound gravitational 
orbits. The particles thereafter become concentrated near the center of such objects, and annihilate more rapidly. 
The neutrinos from the annihilation products can be searched for in neutrino telescopes 

Neutralino annihilations in the galactic halo proceed at a very small rate, but the annihilation products are very 
energetic and might thus be detectable. These are (focusing on antiparticles) gamma rays, antiprotons, positrons, 
and antideuterons. Monochromatic gamma rays are particularly interesting; the \X ~^ 77 > -^°7 branching ratios are 
typically quite small, but the signal is unmistakable and has no astrophysical background. 
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Figure 4: Direct detection of neutralinos. The spin-dependent elastic scattering cross sections of neutralinos on protons and 
neutrons are plotted against neutralino mass. Symbols are as in figure 



6. CAN WE LIVE WITHOUT DARK MATTER? 



Observed matter can not account for the gravitational potentials of objects at the dwarf galaxy scale and larger, thus 
dark matter was hypothesized. An alternative possibility is that gravity weakens less quickly at large distances, or as 
proposed by Milgrom at small acceleration scales. The latter idea is known as MOND, for MOdified Newtonian 
Dynamics. Basically, imposing a minimum acceleration scale an, ~ cHq turns out to explain spiral galaxies quite 
well. If this MOND scale an is allowed to run, clusters might be explained as well. Because gravitational lensing 
mass estimates agree with dynamical estimates, the MOND acceleration scale must apply to photons as well. A 
relativistic theory is needed in any case to construct a consistent cosmology. A number of attempts have been made, 
each with fatal inconsistencies [58j. Quite recently, Bekenstein has proposed a tensor- vector-scalar theory (TeVeS) 
that is successful in reproducing MOND in the proper limit even for photons, and respecting the classical tests of 
general relativity • Whether a consistent cosmology can be constructed with the theory remains to be seen. Note 
that even this is a "dark matter" theory, in that the structure of galaxies is still explained by the introduction of 
new particles. 



7. DISCUSSION 



We are unable to escape the fact that new physics is required to explain the large scale structure of the universe. 
Measurements of the cosmological matter density, most recently by the WMAP satellite, are precise indications (7% 
error, and fast improving) of something relating to this new physics. Identifying exactly what that something is is 
one of the chief challenges in physics in the next decade(s). 
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